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'• Abstract 

> 

O I Using data collected with the LHCb detector in proton-proton collisions at a centre- 

I of- mass energy of 7 TeV, the hadronic decay J/'?/;/o(980) is observed. This 

^ ■ CP eigenstate mode could be used to measure mixing-induced CP violation in the 

O ' system. Using a fit to the tt+tt^ mass spectrum with interfering resonances gives 

Rfo/^ = = 0.252t'o'otlToll In the interval ±90 MeV around 

980 MeV, corresponding to approximately two full /o widths we also find R' = 
^ ; r(B,^j/t^7r^Tr^^^rn^7T^7r^M80^Mcv| <90 Mcv) ^ ^ ^ ^ ^ 0.016, where in both cases 

^ ■ the uncertainties are statistical and systematic, respectively. 
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1 Introduction 



In decays some final states can be reached either by a direct decay amplitude or via 
a mixing amplitude. For the case of J/ '^4' decays, the interference between these 

two amplitudes allows observation of a CP violating phase. In the Standard Model (SM) 
this phase is —2(3s = — 0.036^Q;oQ^g radians, where (3s = aTg{—VtsV^l/VcsV*fj), and the Vij 
are CKM matrix elements [1]. This is about 20 times smaller in magnitude than the 
measured value of the corresponding phase 2/3 in B^ mixing. Being small, this phase can 
be drastically increased by the presence of new particles beyond the SM. Thus, measuring 
Ps is an important probe of new physics. 

Attempts to determine Ps have been made by the CDF and DO experiments at the 
Tevatron using the B^ — )■ J/ipcf) decay mode [2]. While initial results hinted at possible 
large deviations from the SM, recent measurements are more consistent [3lll]. However, 
the Tevatron limits are still not very constraining. Since the final state consists of two 
spin-1 particles, it is not a CP eigenstate. While it is well known that CP violation can be 
measured using angular analyses [S] , this requires more events to gain similar sensitivities 
to those obtained if the decay proceeds via only CP-even or CP-odd channels. In Ref. |6] 
it is argued that in the case of J/ipct) the analysis is complicated by the presence of an 
S-wave K'^K" system interfering with the that must be taken into account, and that 
this S-wave would also manifest itself by the appearance of /o(980) — )■ tt+tt" decays. This 
decay B^ — > J / ip fo{980) is to a single CP-odd eigenstate and does not require an angular 
analysis. Its CP violating phase in the Standard Model is —2/3^ (up to corrections due to 
higher order diagrams). In what follows, we use the notation /o to refer to the /o(980) 
state. 

By comparing D'^ — /ovr"*" decays where the /o was detected in both K^K~ and tt+tt" 
modes it was predicted that [B] 

_ r(i?° ^ J/^/o, /o ^ TT+TT-) 

A decay rate at this level would make these events very useful for measuring [5s if back- 
grounds are not too large. 

The dominant decay diagram for these processes is shown in Fig. [T] It is important 
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Figure 1: Decay diagram for — )■ J/ip{fo or 0) decays, 
to realize that the ss system accompanying the J/ip is an isospin singlet (isoscalar), and 



thus cannot produce a single meson that is anything but isospin zero. Thus, for example, 
in this spectator model production of a p'^ meson is forbidden. The dominant low mass 
isoscalar resonance decaying into tt+tt^ is the /o(980) but other higher mass objects are 
possible. 

Although the /q mass is relatively well estimated at 980±10 MeV (we use units with 
c = 1) by the PDG, the width is poorly known. Its measurement appears to depend on 
the final state, and is complicated by the opening of the KK channel close to the pole; 
the PDG estimates 40—100 MeV [H]. Recently CLEO measured these properties in the 
semileptonic decay — )■ foc^v, where hadronic effects are greatly reduced, determining 
a width of (9ll^^ ± 3) MeV [12]. 

2 Data sample and analysis requirements 

We use a data sample of approximately 33 pb~^ collected with the LHCb detector in 
2010 |7j. The detector elements are placed along the beam line of the LHC starting 
with the Vertex Locator (VELO), a silicon strip device that surrounds the proton-proton 
interaction region and is positioned 8 mm from the beam during collisions. It provides 
precise locations for primary pp interaction vertices, the locations of decays of long-lived 
particles, and contributes to the measurement of track momenta. Other devices used to 
measure track momenta comprise a large area silicon strip detector (TT) located in front 
of a 3.7 Tm dipole magnet, and a combination of silicon strip detectors (IT) and straw 
drift chambers (OT) placed behind. Two Ring Imaging Cherenkov (RICH) detectors are 
used to identify charged hadrons. Further downstream an Electromagnetic Calorimeter 
(ECAL) is used for photon detection and electron identification, followed by a Hadron 
Calorimeter (HCAL), and a system consisting of alternating layers of iron and chambers 
(MWPC and triple-GEM) that distinguishes muons from hadrons (MUON). The ECAL, 
MUON, and HCAL provide the capability of first-level hardware triggering. 

This analysis is restricted to events accepted hy & J/ifj — )■ trigger. Subse- 

quent analysis selection criteria are applied that serve to reject background, yet preserve 
high efficiencies on both the J/t/^vt+vt" and J /%pK^ K~ final states, as determined by 
Monte Carlo events generated using PYTHIA [S], and LHCb detector simulation based 
on GEANT4 [5]. Tracks are reconstructed as described in Ref. [7]. To be considered as 
a, J/%1) ^ fi'^fi~ candidate opposite sign tracks are required to have transverse momen- 
tum, pt, greater than 500 MeV, be identified as muons, and form a common vertex with 
fit psr number of degrees of freedom (ndof) less than 11. The fi~^fi~ invariant mass 
distribution is shown in Fig. |5] with an additional requirement, used only for this plot, 
that the pseudo proper-time, t^, be greater than 0.5 ps, where is the distance that the 
J/ip candidate travels downstream parallel to the beam, along z, times the known J/ip 
mass divided by the z component of the candidate's momentum. The data are fit with a 
Crystal Ball signal function [10] to account for the radiative tail towards low mass, and a 
linear background function. There are 549,000±1100 J/ip signal events in the entire mass 
range. For subsequent use only candidates within ±48 MeV of the known J/ip mass are 
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selected. 




Figure 2: The invariant mass for candidates satisfying the trigger and analysis 

requirements and having > 0.5 ps. The data points are shown as circles; the error 
bars are smaller than the circle radii. The dashed line shows the Crystal Ball signal 
function [10], the dotted line the background and the solid line the sum. 

Pion and kaon candidates are selected if they are inconsistent with having been pro- 
duced at the closest primary vertex. The impact parameter (IP) is the minimum distance 
of approach of the track with respect to the primary vertex. We require that the formed 
by using the hypothesis that the IP is equal to zero be > 9 for each track. For further con- 
sideration these tracks must be positively identified in the RICH system. Particles forming 
opposite-sign di-pion candidates must have their scalar sum > 900 MeV, while those 
forming opposite-sign di-kaon candidates must have their vector sum > 1000 MeV, 
and have an invariant mass within ±20 MeV of the mass. 

To select candidates we further require that the two pions or kaons form a vertex 
with a < 10, that they form a candidate vertex with the J/ip where the vertex fit 
X^/ndof < 5, and that this B^ candidate points to the primary vertex at an angle not 
different from its momentum direction by more than 0.68°. 

Simulations are used to evaluate our detection efficiencies. For the J/ ipcf) final state we 
use the measured decay parameters from CDF [3] . The J / ipfo final state is simulated using 
full longitudinal polarization of the J/ ip meson. The efficiencies of having all four decay 
tracks in the geometric acceptance and satisfying the trigger, track reconstruction and 
data selection requirements are (1.471±0.024)% for J/ipfo, requiring the ti^ti~ invariant 
mass be within ±500 MeV of 980 MeV, and (1.454±0.021)% for J/^(p, having the K+K' 
invariant mass be within ±20 MeV of the mass. The uncertainties on the efficiency 
estimates are statistical only. 
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3 Results 



The J/ipK~^K~ invariant mass distribution is shown in Fig. |3l The di-muon invariant 
mass has been constrained to have the known value of the J/ip mass; this is done for 
all subsequent invariant mass distributions. The data are fit with a Gaussian signal 
function and a linear background function. The fit gives a mass of 5366. 7±0. 4 MeV, 
a width of 7.4 MeV r.m.s., and a yield of 635±26 events. 
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Figure 3: The invariant mass of J /ipK^K' combinations when the K^K~ pair is required 
to be with ±20 MeV of the mass. The data have been fit with a Gaussian signal function 
whose mass and width are allowed to float and linear background function shown as a 
dashed line. The solid curve shows the sum. 

Initially, to search for a /o(980) signal we restrict ourselves to an interval of ±90 MeV 
around the /o mass, approximately two full /o widths [12]. The candidate invariant 
mass distribution for selected J/ipn~^7r~ combinations is shown in Fig. m The signal 
is fit with a Gaussian whose mean and width are allowed to float. We also include a 
background component due to B^ — ?■ J/'ipn^n" that is taken to be Gaussian, with mass 
allowed to float in the fit, but whose width is constrained to be the same as the B^ signal. 
Other components in the fit are B^ — )■ J/i/jK*^, combinatorial background taken to have 
an exponential shape, B~^ — )■ J /%l)K^{oi 7r+), and other specific decay backgrounds 
including 5^ — )■ J/iprj', 1]' — )■ p7, B^ — )■ J/'ip4', (f) ~^ tt^tt'tt^. The shape of the sum of 
the combinatorial and B^ — t- J /ipK^iji'^) components is taken from the like-sign events. 
The shapes of the other components are taken from Monte Carlo simulation with their 
normalizations allowed to float. 

We perform a simultaneous unbinned likelihood fit to the tt'^tt" opposite-sign and sum 
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of 7r"'"7r+ and tt tt like-sign event distributions. The fit gives a mass of 5366.1±1.1 
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Figure 4: (a) The invariant mass of J/ tpTT^ir" combinations when the tt+tt" pair is required 
to be with ±90 MeV of the /o(980) mass. The data have been fit with a signal Gaussian 
and several background functions. The thin (red) solid curve shows the signal, the long- 
dashed (brown) curve the combinatorial background, the dashed (green) curve the — )■ 
J/xl:K~^{n~^) background, the dotted (blue) curve the — J/ipK*^ background, the 
dash-dot curve (purple) the B^ J / ipTi'^ tx~ background, the barely visible dotted curve 
(black) the sum of 5° — > J / iprj' and J/ ipcj) backgrounds, and the thick-solid (black) curve 
the total, (b) The same as above but for like-sign di-pion combinations. 



MeV in good agreement with the known mass of 5366. 3±0. 6 MeV, a Gaussian width 
of 8.2±1.1 MeV, consistent with the expected mass resolution and 111±14 signal events 
within ±30 MeV of the 5° mass. The change in twice the natural logarithm of the 
fit likelihood when removing the -B^ signal component, shows that the signal has an 
equivalent of 12.8 standard deviations of significance. The like-sign di-pion yield correctly 
describes the shape and level of the background below the -B° signal peak, both in data 
and Monte Carlo simulations. There are also 23±9 B^ — t- J/iI)'k'^ti~ events. 
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Having established a clear signal, we perform certain checks to ascertain if the struc- 
ture peaking near 980 MeV is a spin-0 object. Since the 5° is spinless, when it decays 
into a spin-1 J/ip and a spin-0 /o, the decay angle of the J/ip should be distributed as 
1 — cos^Oj/^,, where 6j/^ is the angle of the /i"*" in the J/ip rest frame with respect to 
the direction. The polarization angle, Of^, the angle of the vr"*" in the /o rest frame 
with respect to the direction, should be uniformly distributed. A simulation of the 
J/ip detection efficiency in these decays shows that it is approximately independent of 
cosOj/^. The acceptance for /o — t- tt+tt^ as a function of the vr"^ decay angle shows an 
inefficiency of about 50% at cosOf^^ = ±1 with respect to its value at cosOf^ = 0. It is fit 
to a parabola and the inefficiency corrected in what follows. 

The like-sign background subtracted J/ip helicity distribution is fit to a 1 — acos^ 9j/^ 
function as shown in Fig. [5]^a). The fit gives a = 0.81 ± 0.21 consistent with a longitu- 
dinally polarized J/ip (spin perpendicular to its momentum) and a spin-0 /o meson. The 

of the fit is 10.3 for 8 degrees of freedom. Similarly, we subtract the like-sign back- 




Figure 5: Angular distributions of events within ±30 MeV of the B^ mass and ±90 MeV 
of the /o mass after like-sign background subtraction, (a) The cosine of the angle of the 
/i"*" with respect to the B^ direction in the J/ip rest frame for B^ — j- J/?/^7r"*"7r~ decays. 
The data are fit with a function f{cos6j/^) = 1 — acos^ (^j/^p- (b) The cosine of the angle 
of the 7r~^ with respect to the B^ direction in the di-pion rest frame for B^ — )■ J/t/^tt+tt" 
decays. The data are fit with a fiat line. 

ground and fit the efficiency corrected vr"''-/?" helicity distribution to a constant function as 
shown in Fig. |5]^b). The fit has a x^/ndof equal to 15.9/9, still consistent with a uniform 
distribution as expected for a spinless particle. 

To view the spectrum of tt'^tt' masses, between 580 and 1480 MeV, in the J/t/^vt+vt" 
final state we select events within ±30 MeV of the B^ and plot the invariant mass spectrum 
in Fig. O The data show a strong peak near 980 MeV and an excess of events above the 
like-sign background extending up to 1500 MeV. Our mass spectrum is similar in shape to 
those seen previously in studies of the S-wave tt+tt" system with s's quarks in the initial 
state p^] - [T5] . To establish a value for -R/o/0 requires fitting the shape of the /o resonance. 
Simulation shows that our acceptance is independent of the 7i~^7i~ mass, and we choose 
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Figure 6: The invariant mass of tt+tt^ combinations when the J j'i\}TX^TX~ is required to be 
within ±30 MeV of the i?^ mass. The dashed curve is the hke-sign background that is 
taken from the data both in shape and absolute normahzation. The dotted curve is the 
result of the fit using Eq. H] and the solid curve the total. 



an interval between 580 and 1480 MeV. Guidance is given by the BES collaboration who 
fit the spectrum in J/-?/; — )■ 07r+7r~ decays [2]. We include here the /o(980) and /o(1370) 
resonances, though other final states may be present, for example the /2(1270) a 2"*""'" 
state [13l[Tl] ; it will take much larger statistics to sort out the higher mass states. We use 
a coupled-channel Breit-Wigner amplitude (Platte) for the /o(980) resonance [16] and a 
Breit-Wigner shape (BW) for the higher mass /o(1370). Defining m as the tt+tt" invariant 
mass, the mass distribution is fit with a function involving the square of the interfering 
amplitudes 

\A{m)f = Nomp{m)q{m) |Flatte[/o(980)] + exp^'^) BW[/o(1370)] f , (2) 

where A'^o is a normalization constant, p{m) is the momentum of the 7r+, q{m) the mo- 
mentum of the J/ip in the tt+tt" rest-frame, and 6 is the relative phase between the two 
components. The Platte amplitude is defined as 

Platte(m) = — ■, — ] ■ , (3) 

mg - - imo{giPnn + g2PKK) 

where mo refers to the mass of the /o(980) and Pt^tt and pkk sue Lorentz invariant phase 
space factors equal to 2p{m)/m for p^^r- The g2PKK term accounts for the opening of the 
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kaon threshold. Here pkk = 2px{m)/m where pxim) is the momentum a kaon would 
have in the Tr"*"?!" rest-frame. It is taken as an imaginary number when m is less than 
twice the kaon mass. We use mogi = 0.165 ± 0.018 GeV^, and 5'2/fi'i = 4.21 ± 0.33 as 
determined by BES [H]. 

The /o(1370) mass and width values used here are 1434±20 MeV, and 172±33 MeV 
from an analysis by E791 [13]. We fix the central values of these masses and widths 
in the fit, as well as mogi and the 5'2/fi'i ratio for the /o(980) amplitude. The mass 
resolution is incorporated as a Gaussian convolution in the fit as a function of tt'^tt" 
mass. It has an r.m.s. of 5.4 MeV at 980 MeV. We fit both the opposite-sign and like-sign 
distributions simultaneously. The results of the fit are shown in Fig. [HI The x^/ndof 
is 44/56. We find an /o(980) mass value of 972±25 MeV. There are 265±26 events 
above background in the extended mass region, of which (64l^g)% are associated with 
the /o(980), (12 ±4)% are ascribed to the /o(1370) and (24^1)% are from interference. 
The fit determines 5 = 61 ± 36°. The fit fraction is defined as the integral of a single 
component divided by the coherent sum of all components. The /o(980) yield is 169^21 
events. The lower mass cutoff of the fit region loses 1% of the /o(980) events. The change 
in twice the log likelihood of the fit when removing the /o(980) component shows that it 
has an equivalent of 12.5 standard deviations of significance. 

Using the 169 /o events from J / ipTi^ ti~ , and the 635 (f) events from J/'iIjK~^K~ , cor- 
recting by the relative efficiency, and ignoring a possible small S-wave contribution under 
the peak [T7j, yields 

D ^ r(E° ^ J/^/o, /o ^ TT+TT-) _ +0.046+0.027 

- r(50 ^ ^ - U-^52_o.o32-0.033 • W 

Here and throughout this Letter whenever two uncertainties are quoted the first is statisti- 
cal and the second is systematic. This value of depends on the decay amplitudes used 
to fit the 7r"'"7r~ mass distribution and could change with different assumptions. To check 
the robustness of this result, an incoherent phase space background is added to the above 
fit function. The number of signal /o(980) events is decreased by 7.3%. If we leave the 
/o(1370) out of this fit, the original /o(980) yield is decreased by 6.5%. The larger number 
of these two numbers is included in the systematic uncertainty. The BES collaboration 
also included a a resonance in their fit to the vr+vr" mass spectrum in J/ip — )■ ^tt+vt" 
decays [2]. We do not find it necessary to add this component to the fit. 

The systematic uncertainty has several contributions listed in Tabled! There is an un- 
certainty due to our kaon and pion identification. The identification efficiency is measured 
with respect to the Monte Carlo simulation using samples of D*~^ — )■ tt^D^, — K~7r~^ 
events for kaons, and samples of Kg — )■ tt+tt" decays for pions. The correction to Rf^/^ is 
0.947±0.009. This correction is already included in the efficiencies quoted above, and the 
1% systematic uncertainty is assigned for the relative particle identification efficiencies. 

The efficiency for detecting </> — > K^K~ versus a tt+tt" pair is measured using 
meson decays into ^vr"*" and J^'~7r+7r+ in a sample of semileptonic B decays where 
B D^Xn~V [18]. The simulation underestimates the efficiency relative to the tt+vt" 
efficiency by (6±9)%, so we take 9% as the systematic error. 
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Besides the sources of uncertainty discussed above, there is a variation due to varying 
the parameters of the two resonant contributions. We also include an uncertainty for a 
mass dependent efficiency as a function of vr+vr" mass by changing the acceptance function 
from flat to linear and found that the /o yield changed by 2.3%. The difference AF/F 
between CP even and CP odd Bg eigenstates is taken as 0.088. Ignoring this difference 
results in less than a 1% change in the relative efficiency. 



Table 1: Relative systematic uncertainties on -R/o/t^ (%)• Both negative and positive 
changes resulting from the parameter variations are indicated in separate columns. 



Parameter 


Negative change 


Positive change 


/o(1370) mass 


0.3 


1.9 


/o(1370) width 


2.3 


2.6 


7r+7r~ mass dependent efficiency 


2.3 


2.3 




4.2 


3.6 


92/91 


0.7 


0.7 


Addition of non-resonant tt+tt^ 


7.3 





MC statistics (efficiency ratio) 


2.3 


2.3 


5° Pt distribution 


0.5 


0.5 


5° mass resolution 


0.5 


0.5 


PID efficiency 


1.0 


1.0 


(j) detection 


9.0 


9.0 


Total 


13.1 


10.8 



In order to give a model independent result we also quote the fraction, i?', in the 
interval ib90 MeV around 980 MeV, corresponding to approximately two full-widths, 
where there are 111±14 events. Then 

H . r(B;^.W.^.-. |M.^.-)-980 Mey|<90MeV) _ 

F(50 ^ J/^(t), ^ K+K-) 

(5) 

This ratio is based on the fit to the 5° mass distribution and does not have any uncer- 
tainties related to the fit of the tt+tt" mass distribution. Based on our fits to the tt+tt" 
mass distribution, there are negligible contributions from any other signal components 
than the /o(980) in this interval. 

The original estimate from Stone and Zhang was -R/o/t^ = 0.20 [6]. More recent pre- 
dictions have been summarized by Stone [19] and have a rather wide range from 0.07 to 
0.50. 



4 Conclusions 

Based on the polarization and rate estimates described above, the first observation of a 
new CP-odd eigenstate decay mode of the meson into J /ipf o{980) has been made. 
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Using a fit including two interfering resonances, the /o(980) and /o(1370), the ratio to 
J /iljcj) production is measured as 

p = ^ J/i^h. /o ^TT+TT-) _ +0.046+0.027 ( a^-s 

Hf,/^ - ^^^0 ^ J/^^^ ^ ^ K+K~) ~ ^-^^^-0.032-0.033 ■ Ibj 

By selecting events within ±90 MeV of the /o(980) mass the ratio becomes R' = 0.162 ± 
0.022 ±0.016 . 

The events around the /o(980) mass are large enough in rate and have small enough 
backgrounds that they could be used to measure (3s without angular analysis. It may also 
be possible to use other data in the tt+tt^ mass region above the /o(980) for this purpose 
if they turn out to be dominated by S-wave. 
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